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Abstract 

A new species of the Stegana (Steganina) ornatipes species group (Diptera: Drosophilidae) is 
described from Hainan, China, S. (S.) xipengi sp. nov. Based on the mitochondrial ND2 and COI 
gene sequences, the relationships among eight species from mainland China of the ornatipes 
group, and their relationships to the undulata, nigrolimbata and shirozui species groups of the 
same subgenus, are investigated, using two species of the subgenus Stegana, S. emeiensis and S. 
quadrata, as outgroups. The result shows that S. (S.) mengla is debarred from the ornatipes 
group. 
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So far five species groups have been identified 
in the subgenus Steganina Wheeler (Diptera: 
Drosophilidae) of the genus Stegana Meigen: 
coleoptrata group (Lastovka and Maca 1982; 
Chen and Chen 2008), nigrolimbata group 
(Sidorenko 2002; Cao and Chen 2008), 
shirozui group (Chen et al. 2009), undulata 
group (Sidorenko 2002) and ornatipes group 
(Cheng et al. 2009), and they included 51 
species; most of them were from the Oriental 
region except for some species of the 
coleoptrata group from the Palearctic region. 
The ornatipes group includes ten species from 
the Oriental region: S. (S.) vietnamensis 
Sidorenko, 1997 from Virtnam; S. (S.) 
albiventralis Cheng, Gao et Chen, 2009; S. 
(S.) angusigena Cheng, Gao et Chen, 2009; S. 
(S.) chitouensis Sidorenko, 1998; S. (S.) 
lingnanensis Cheng, Gao et Chen, 2009; S. 
(S.) mengla Cheng, Gao et Chen, 2009; S. (S.) 
nulliseta Cheng, Gao et Chen, 2009; S. (S.) 
ornatipes Wheeler et Takada, 1964; S. (S.) 
pilosella Cheng, Gao et Chen, 2009 and 5*. (S.) 
zhaofengi Cheng, Gao et Chen, 2009 from 
China. This group is supported by the 
following morphological characters as the 
diagnosis: surstylus large, with a strong 
prensiseta apically and several thin, long 
setae; 10th sternite mostly narrowed, nearly 
arcuate, with a pair of projections 
posterolaterally; gonopods with a pair of 
projections sublaterally. On the other hand, 
this group is similar to the nigrolimbata group 
in sharing the following morphological 
characters: palpus mostly black, sometimes 
yellow basally; gena yellow to brown, narrow 
(ch/o < 0.10); aedeagus basally contiguous to 
aedeagal apodeme; which shows the both are 
more closely related to each other than other 
members of the subgenus Steganina. 
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Recently, some studies of molecular 
phylogeny were appeared to the subfamily 
Steganinae (Otranto et al. 2008; He et al. 
2009a, b; Zhao et al. 2009; Li et al. 2010). 
Otranto et al. (2008) reconstructed the 
phylogenetic relationships among 13 species 
of 8 genera of Steganinae based on the DNA 
sequences of the cytochrome oxidase subunit I 
{COT) gene, however, in their phylogenetic 
analysis, only two Stegana species were 
sampled as the representative. Li et al. (2010) 
investigated the phylogenetic relationships 
among seven of the Chinese species of the 
subgenus Stegana (s.s.) based on the DNA 
sequences of the NADH dehydrogenase 
subunit 2 (ND2) gene, using two species of 
the subgenus Steganina (S. nigrilimbata Duda, 
1924 and S. ctenaria Nishiharu, 1979) as 
outgroup taxa. 

In the present study, we described a new 
species of the ornatipes group from Hainan, 
China. We also constructed the molecular 
phylogeny based on the mtDNA sequences of 
ND2 and COI genes. To investigate the 
relationships in ornatipes group and with the 
other species groups of subgenus Steganina, 
we employed the additional seven species 
from mainland China of this group and S. 
nigrolimbata Duda, S. xiaoleiae Cao and 
Chen, S. ctenaria Nishiharu, and S. undulata 
de Meijere which belong to nigrolimbata, 
shirozui and undulata species groups of the 
subgenus Steganina as ingroup taxa Two 
species from subgenus Stegana, S. emeiensis 
Sidorenko and S. quadrata Cao and Chen 
were chosen as outgroup taxa. 

Materials and Methods 

All materials were collected on tussock and 
tree trunks along streams in forest, preserved 
in 75% ethanol immediately and identified 
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(Table 1). A small piece of tissue was 
removed from the fly abdomen and used for 
the DNA extraction; then, the body and 
terminalia parts were dried and deposited in 
the Department of Entomology, South China 
Agricultural University, Guangzhou, China 
(SCAU). McAlpine (1981) was followed for 
morphological terminology and Zhang and 
Toda (1992), and Chen and Toda (2001) for 
the definitions of measurements, indices and 
abbreviations. 

DNA extraction and sequencing 

The total DNA was extracted using the DNA 
extraction Kit (TIANGEN®) according to the 
manufacture's protocol. The ND2 gene and 
the 5' end of COI gene were amplified. 
Primers used were given in table 2. The PCR 
cycle program comprised an initial 3 min of 
predenaturation at 94 °C, 35 cycles of 
amplification ( 50 s of denaturation at 94 °C; 1 
min of annealing at 53 °C for ND2, 49 °C for 
COI; 1 min of extension at 72 °C), and a final 
elongation for 5 min at 72 °C. When possible, 
purified amplified products were directly run 
on an ABI 3730 sequencer for sequenceing, 
otherwise they were cloned into the pMD18-T 
plasimid vector (TAKARA®), and then 
sequenced. The related ND2 sequences of S. 
emeiensis, S. quadrata, S. ctenaria and S. 
nigrolimabata were retrieved from the 



National Center for Biotechnology 
Information (NCBI); the related COI 
sequences of emeiensis, S. nigrolimbata and 5*. 
undulata were also retrieved from the NCBI. 

Phylogenetic analyses 

The sequences were aligned by the Clustal W 
(Thompson et al. 1994) method implemented 
in program MEGA 4.0 (Tamura et al. 2007) 
with default options. A partition homogeneity 
test (PHT) between the ND2 and COI 
sequences was performed with PAUP 4.0b 10* 
(Swofford 2002). The program DAMBE 
5.0.80 (Xia and Xie 2001) was used to 
measure the nucleotide substitution saturation 
using the method of Xia et al. (2003) as the 
substitution saturation masked the 
phylogenetic signal (Lopez et al. 1999; 
Philippe and Froterre 1999). Base 
compositions were investigated by means of 
the software PAUP 4.0b 10* (Swofford 2002), 
and a x 2 test was also used to test the 
nucleotide composition homogeneity. 
Uncorrected pairwise divergence was 
estimated by program MEGA 4.0 (Tamura et 
al. 2007). 

Phylogenetic trees were constructed by using 
the maximum parsimony (MP) and maximum 
likelihood (ML) in PAUP 4.0b 10* (Swofford 
2002), the Bayesian inferring (BI) method 



Table I. Collecttion data of samples for DNA sequencing, and accession numbers of the ND2 and COI sequences. 



Subgenera 


Groups 


Species 


Species 
numbers in 
SCAU 


Localities 


Latitude and longitude 


Accession numbers of ND2 


Accession numbers of COI 


Stegana 




emeiensis 


1 20268 


Menglun, Yunnan, China 


2I"4I'N, I0I"25'E 


EU8055 1 5 


HM636455 






quadrata 


1 20084 


Mengla, Yunnan, China 


2I"28'N, 101 "38'E 


EU8055 1 6 


HQ270I47 


Steganina 


nigrilimbata 


nigroiimbata 


1 20390 


Guangzhou, Guangdong, China 


23'10'N, 1 I2"34'E 


EU8055 1 3 


HM636458 






xiaoleiae 


1 20453 


Maoershan, Guangxi, China 


25"5I'N, 1 I0"27'E 


GQ259982 


HQ260633 




shltozul 


ctenaria 


120321 


Kumamoto, Japan 


32"37'N, I30"5I'E 


EU8055I4 


HQ270I48 




undulata 


undulata 


120175 


Mengla, Yunnan, China 


Ditto 


GQ249I96 


HM646466 




ornatipes 


albiventralis 


120610 


Jingdong, Yunnan, China 


24'32'N, 1 0 1 "0 1 E 


GQ259983 


HQ260634 






angustigena- GX 


1 20627 


Nonggang, Guangxi, China 


22'22'N, 106'51'E 


GQ259984 


HQ260637 






angustigena-HN 


1 20626 


Jianfengling, Hainan, China 


18'41'N, I08"52'E 


GQ259985 


HQ260635 






angustigena-YN 


1 20634 


Mengla, Yunnan, China 


Ditto 


GQ259986 


HQ260636 






lingnanensis 


120636 


Guangzhou, Guangdong, China 


Ditto 


GQ259987 


HQ260638 






pilosella 


1 20642 


Mengla, Yunnan, China 


Ditto. 


GQ259988 


HQ260640 






nulliseta 


120656 


Mengla, Yunnan, China 


Ditto 


GQ25999I 


HQ260639 






xipengi sp. nov. 


120691 


Jianfengling, Hainan, China 


Ditto 


GQ259992 


HQ26064I 






zhaofengi-YN 1 


1 20459 


Mengla, Yunnan, China 


Ditto. 


GQ259993 


HQ260643 






z/ioofengi-YN2 


1 20464 


Mengla, Yunnan, China 


Ditto 


GQ259994 


HQ260642 






zf)<jofengi-YN3 


1 20665 


Mengla, Yunnan, China 


Ditto 


GQ259995 


HQ260644 




7 


meng/a-YN 1 


120653 


Mengla, Yunnan, China 


Ditto 


GQ259989 


HQ260645 






meng/o-YN2 


120652 


Mengla, Yunnan, China 


Ditto 


GQ259990 


HQ260646 
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Target gene 


Primer 


Primer sequence (5"-3') 


Reference 


ND2 


ND2-H 


AAGCTACTGGGTTCATACC 


Park 1999 




ND2-T3 


AGGCGATAGATTGTAAATC 


Li 2010 


CO/ 


COI-FI 


CGCCTAAACTTCAGCCACTT 


He et al. 2009a 




LCOI490 


GGTCAACAAATCATAAAGATATTGG 


Folmer et al. 1994 




HC02I98 


TAAACTTCAGGGTGACCAAAAAATCA 


Folmer et al. 1994 



performed in MrBayes 3.2.1 (Huelsenbeck 
and Ronquist 2001; Ronquist and 
Huelsenbeck 2003). The MP and ML trees 
were searched by the heuristic method, with 
initial trees obtained by randomly adding taxa, 
and the TBR algorithm was used in branching 
swapping. Branch support for each node in the 
MP and ML trees was assessed by 1000 
bootstrap replicates. The nucleotide 
substitution models of ML and BI analyses 
were selected by MrModeltest 2.3 (Nylander 
2004) using the hierarchical likelihood ratio 
test (hLRT) criterion (Posada and Crandall 
1998). In the BI analyses, the site-specific 
models were assigned to dataset partitioned by 
locus (2 data partitions) and by codon 
positions (6 data partitions). Two independent 
runs with 2,000,000 generations were 
implemented in parallel, sampling frequency 
of every 100 generations was employed. 
When the average deviation of split 
frequencies fell well below 0.01, the two runs 
were stopped. For each running, the 5,000 
early-phase samples were burn-in, the rest 
samples were used in summarizing and a 
majority rule tree showing all the compatible 
partitions was obtained. 



lends conflict to a given node, and zero 
indicate that the partition lends neither support 
nor conflict to a given node. The partitioned 
Bremer support values were calculated using 
the partitioned constraint file in TreeRot v3 
(Sorenson 1999). 

Results 

Stegana (Steganina) xipengi sp. nov. 
(Figures 1 , 2) 

Diagnosis 

This species is related to S. (S.) albiventralis 
from Yunnan in having the entirety white 
katepisternum, but clearly distinguishable 
from it by the palpus yellow basally, black 
distally, the mesonotum brown, without stripe 
(in albiventralis: palpus entirely yellow; 
mesonotum brown, with yellow stripe 
medially). 

Description 

Male: Frons and face not rectangular in 
profile. Eyes red. Ocellar triangle black, with 
1 pair of small setae above ocellar setae. 



Partition Bremer support (PBS) was used to 
show the contribution of each gene partition to 
the Bremer support of the simultaneous 
analysis (Baker and DeSalle 1997). Values 
can be positive, negative or zero and sum of 
all the partitioned Bremer support values at a 
node will equal the Bremer support value for 
that node. A positive PBS value suggests 
support for the node by that gene, whereas a 
negative PBS value indicates that the partition 




Figures I -2. Stegana (Steganina) xipengi sp. nov., S- I. 
Epandrium, cercus and surstylus (lateral view); 2. surstylus 
(frontal view). Scale bars = 0. 1 mm. High quality figures are 
available online. 
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Postvertical setae slightly behind vertex ridge. 
Frons shiny, brown, with sporadic, minute 
setulae submedially, and a black, transverse 
band above ptilinal fissure. Proclinate orbital 
setae slightly nearer to ptilinal fissure than to 
inner vertical setae. Pedicel brown; first 
flagellomere yellow only basally, mostly 
black. Face black with yellow, transverse 
band medially, broadened ventrally; facial 
carina absent. Clypeus black medially, yellow 
laterally. Palpus yellow basally, black distally, 
with 1-2 longer setae distally and several 
shorter setae basally. Gena yellow, narrow. 
Vibrissa prominent; other orals small. Occiput 
glossy, yellow, but black around occipital 
foramen. Mesonotum brown. Mesopleuron 
with a black longitudinal stripe above 
(running from propleuron to base of halter). 
Postpronotal lobe brown on upper part, white 
on lower part, with 1 long and a few short 
setae. Acrostichal setulae approximately in 10 
irregular rows. Prescutellar setae 1 pair. 
Katepisternum entirely white. Scutellum 
brown; basal setae divergent; apical setae 
crossing with each other. Wing dark brown 
anteriorly, pale posteriorly, curved downward 
on distal part. Basal medial-cubital crossvein 
present. Ci with 2 isometric setae. Costal vein 
with 9 minute spinules on ventral surface 
between veins R2+3 and Pm +5 . Vein R2+3 
obviously curved to costa at tip; Veins R4+5 
and Mi convergent distally. Halters white 
basally, greyish brown distally. Legs whitish 
yellow, brown on apical part of fore femur, 
and fore and hind tarsomeres, dark brown to 
black on medially on mid and hind femora, 
with 2 dark brown rings on fore and mid 
tibiae. Fore femur with 3-4 setae on distal part 
of ventral surface. Apical seta present on mid 
tibia. Preapical dorsal setae present on all 
tibiae. Mid tibia (misused to mid tarsus in Cao 
and Chen 2008; Cheng et al. 2009) with 5 
strong setae on basal part of dorsal surface. 
Mid and hind tarsomeres with 2 and 1 row(s) 
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of minute cuneiform setulae on ventral 
surface, respectively; fore and hind 1st 
tarsomeres slightly shorter than the rest 
combined; mid 1st tarsomere longer than the 
rest combined. Abdominal all tergites dark 
brown. Sternites brown; 3rd to 5th broadened; 
6th covered with 5th. Epandrium pubescent 
except for anteroventral margins, with 
approximately 21 setae near posterior margin 
on each side (Figure 1). Cercus separated 
from epandrium, setigerous, lacking 
pubescence (Figure 1). Surstylus separated 
from epandrium, with several thin, long setae 
on inner margin and surface (Figure 2), 
apically strongly curved and with 1 strong 
prensiseta (Figure 2). The hypandrium, 
gonopods, aedeagus and aedeagal apodeme 
were lost when clearing them in KOH 
solution. 

Measurements 

BL = 2.76 mm in holotype; ThL = 1.32 mm; 
WL = 2.58 mm; WW =1.12 mm. Indices: arb 
= 8/7, avd = 0.83, adf = 1.20, flw = 1.80, 
FW/HW = 0.36, ch/o = 0.08, prorb = 1.16, 
rcorb = 0.82, vb = 0.30, del = 0.40, presctl = 
0.60, sctl = 1.80, sterno = 0.90, orbito = 2.20, 
dep = 0.20, sctlp = 1.00, C = 1.86, 4c = 1.22, 
4v = 1.74, 5x = 1.40, ac = 9.33, M = 0.61, 
C3F = 0.66. 

Type 

Holotype: S (SCAU, No. 120589), CHINA: 
Jianfeng, Ledong, Hainan, 18°41'N, 
108°52'E, alt. 750 m, 14.iv.2008, ex tussock, 
X.P. Chen. 

Etymology 

Patronym of the collector Xipeng Chen 
(SCAU). 

Distribution 

China (Hainan). 
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Table 3. Results of model selection, composition homogeneity test and test of substitution saturation. 





ND2 




CO/ 




Ist-CP* 


2nd-CP 


3rd-CP 


All 




Ist-CP 


2nd-CP 


3rd-CP 


All 


Model selection 




















Model selected 


GTR+G 


HKY+G 


HKY+G 


GTR+I+G 




SYM+G 


F8I 


HKY+I+G 


GTR+I+G 


Proportion of invariable sites (1) 


0 


0 


0 


0.297 




0 


0 


0.1 13 


0.6104 


Gamma distribution shape parameter 


0.2678 


0.1055 


0.5701 


0.4971 




0.1581 


- 


0.534 


1.359 


Composition homogeneity test 




















df 


54 


54 


54 


54 




54 


54 


54 


54 


.' 2 


13.00 


2.67 


44.86 


20.45 




3.67 


0.32 


31.94 


7.36 


P 


1.0000 


1.0000 


0.8075 


1.0000 




1.0000 


1.0000 


0.9927 


1.0000 


Test of substitution saturation 




















Iss 


0.1484 


0.0783 


0.3847 


0.1742 




0.0585 


0.0093 


0.4286 


0.3728 


Iss.c (for an asymmetrical tree) 


0.6935 


0.6935 


0.6935 


0.7605 




0.6793 


0.6793 


0.6793 


0.7347 


P 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


0.0000 


0.0000 


Iss.c (for an extreme asymmetrical tree) 


0.4575 


0.4575 


0.4575 


0.5479 




0.4390 


0.4390 


0.4390 


0.5123 


P 


0.0000 


0.0000 


0.003 1 


0.0000 




0.0000 


0.0000 


0.7226 


0.0000 



* CP = codon position; 



Molecular analysis 
Data set analysis 

The alignment was 1739 nucleotide positions 
(1029 for ND2 and 710 for COI, respectively) 
in length. There were end gaps in the ND2 
sequence of S. undulata (sites 1-22) and in the 
COI sequence of S. xiaoleiae (sites 1-33). The 
base composition of ND2 and COI were 
generally AT rich with a mean of 83% and 
69%, respectively. It contained high AT 
contents in the 3rd (94.2% and 94.2%, 
respectively) codon positions. Performance of 
the Chi-square test was showed in table 3. It 
yielded a homogeneous base composition in 
the 7VD2-alignments and COJ-alignments or in 
the separate condon positions of the two 
mitochondrial genes. 

The test of substitution saturation showed that 
the observed index of substitution saturation 
(Iss) for M)2-alignments or for COI- 
alignments was significantly lower than the 
corresponding critical index substitution 
saturation (Iss.c), indicating that there was 
little saturation in our sequences. However, 
when considering partitions separated by 
codon, we identified substitution saturation in 
the third codon position of the COI- 



alignments [Iss = 0.4286 < Iss.c = 0.4390 (for 
an extreme asymmetrical tree, p = 0.72)] 
(Table 3). Since none of the resulted trees of 
the present study are extremely asymmetric, 
there should be little substitution saturation in 
our sequence. 

Table 4 shows the uncorrected pairwise p- 
distances for the ND2 and COI sequences. 
The genetic divergence of ND2 sequences of 
species within the ornatipes group ranged 
from 5.28% to 14.24%, and genetic 
divergence of COI sequences ranged from 
4.28% to 10.49%, however, when we took no 
account of the S. mengla, the upper limits 
would declined to 8.37% and 8.42% for ND2 
and COI, respectively. Within the ornatipes 
group, divergences between S. mengla and 
other species ranged from 9.36% to 15.64% 
for ND2 and from 8.57% to 10.93% for COI, 
whereas the genetic variance among groups 
ranged from 11.65% to 14.34% for ND2 and 
from 8.57% to 1 1.23% for COI. 

Phylogenetic analysis 

The PHT resulted in a p value of 0.062, 
indicating that no significant incongruence 
was found between the ND2 and COI data 
sets. The best-fit models selected for the ML 
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Table 4. Uncorrected pairwise p-distance among the ND2 and COI sequences of the ornatipes species group 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 1 


12 


13 


14 


IS 


16 


17 


18 


19 


1 


S. emeiensis 




0.0798 


0.1241 


0.1270 


0.1064 


0.1 196 


0.1 167 


0.1 182 


0.1 182 


0.1 182 


0.1 152 


0.1 137 


0.1329 


0.1 137 


0.1 137 


0.1 137 


0.1 152 


0.1241 


0.1241 


2 


S. quadrota 


0.0747 




0.1 108 


0.1 137 


0.1064 


0.1 196 


0.121 1 


0.1 270 


0.1241 


0.1241 


0.1196 


0.1226 


0.1241 


0.121 1 


0.1285 


0.1285 


0.1329 


0.1078 


0.1078 


3 


S. nigrolimbola 


0.1464 


0. 1 365 




0.0768 


0.1064 


0.093 1 


0.0857 


0.0916 


0.0886 


0.0886 


0.0857 


0.0916 


0.0871 


0.0916 


0.0916 


0.0916 


0.0945 


0.1108 


0.1 108 


4 


S. xiookiae 


0.1554 


0.1404 


0.0697 




0.1123 


0.1 108 


0.1004 


0.0975 


0.0975 


0.0975 


0.0931 


0.0975 


0.0901 


0.0901 


0.0990 


0.0990 


0.1019 


0.1064 


0.1064 


5 


S. aenoria 


0.1434 


0.1255 


0.1245 


0.1345 




0.0960 


0.1064 


0.1049 


0.1019 


0.1019 


0.1004 


0.1093 


0.1078 


0.0975 


0.1064 


0.1064 


0.1078 


0.1049 


0.1049 


6 


S undulata 


0.1663 


0.1514 


0.1404 


0.1434 


0.1 165 




0.0975 


0.1004 


0.1004 


0.1004 


0.0916 


0.0990 


0.097S 


0.0931 


0.1019 


0.1019 


0.1049 


0.1093 


0.1093 


7 


S. atbiventralis 


0.1534 


0.1375 


0.1076 


0.1106 


0.1255 


0.1534 




0.0458 


0.0428 


0.0428 


0.7390 


0.0591 


0.0561 


0.0635 


0.0576 


0.0576 


0.0606 


0.1049 


0.1049 


8 


S. angustigena-HN 


0. 1 564 


0. 1 384 


0.1006 


0.0996 


0.1325 


0.1414 


0.0677 




0.0044 


0.0044 


0.0665 


0.0620 


0.0502 


0.0473 


0.0650 


0.0650 


0.0694 


0.1049 


0. 1 049 


9 


S. angustigena-YN 


0.1573 


0 1 394 


0.1056 


0.1006 


0. 1 345 


0. 1 424 


00677 


0 0120 




0.0000 


00694 


0.0620 


0.0532 


0.0S02 


0.0635 


0.0635 


0.0679 


0.1034 


0 1034 


10 


S. ongustigena-GX 


0.1584 


0.1404 


0.1066 


0.1016 


0.1355 


0.1434 


0.0687 


0.0129 


0.0010 




0.0694 


0.0620 


0.0532 


0.0502 


0.0635 


0.0635 


0.0679 


0.1034 


0.1034 


II 


S. lingnonensis 


0.1484 


0.1295 


0.1046 


0.1036 


0.1235 


0.1504 


0.0697 


0.0568 


0.0608 


0.0618 




0.0827 


0.0753 


0.0606 


0.0857 


0.0857 


0.0871 


0.0945 


0.0945 


12 


S. nuliseta 


0.1514 


0.1315 


0.1096 


0.1096 


0.1335 


0.1544 


0.0837 


0.0767 


0.0807 


0.0817 


0.0777 




0.0724 


0.0709 


0.0620 


0.0620 


0.0679 


0.0975 


0.0975 


13 


S. pilosella 


0 1534 


0 1394 


0 1036 


0.1 125 


0 1315 


0 1544 


0 0807 


0 0573 


0 0598 


0 0608 


0.0667 


0 0807 




0.0591 


0 0679 


0 0679 


0 0739 


0 1049 


0 1049 


14 


S. xipengi 


0.1554 


0.1325 


0.1026 


0.1076 


0.1265 


0.1504 


0.0677 


0.0767 


0.0548 


0.0538 


0.0588 


0.0817 


0.0598 




0.0798 


0.0798 


0.0842 


0.1034 


0.1034 


15 


S. zhoofengi-YN 1 


0.1414 


0.1285 


0.0946 


0.1066 


0.1245 


0. 1 544 


0.0737 


0.0687 


0.0697 


0.0707 


0.0647 


0.0677 


0.0717 


0.0697 




0.0000 


0.0059 


0.1049 


0.1034 


16 


S. zhaofengr-YN2 


0.1444 


0.1295 


0.0936 


0.1076 


0.1255 


0.1564 


00747 


0.0528 


0.0707 


0.0717 


00657 


0.0687 


0.0737 


0.0707 


0.0050 




0.0059 


01049 


0.1034 


17 


S.z/toofengi-YN3 


0.1404 


0.1275 


0.0936 


0.1056 


0.1235 


0.1544 


0.0727 


0.0667 


0.0687 


0.0697 


0.0637 


0.0667 


0.0707 


0.0687 


0.0040 


0.0010 




0.1049 


0.1034 


18 


S. meng/a-YN 1 


0.1454 


0.1245 


0.1355 


0 1365 


0.1 175 


0.1375 


0.1404 


0.1384 


0.1384 


0.1394 


0.1375 


0.1434 


0.1424 


0.1315 


0.1345 


0.1355 


0.1345 




0.0015 


19 


S. meng/o-YN2 


0.1454 


0.1245 


0.1355 


0.1365 


0.1175 


0.1375 


0.1404 


0.1384 


0.1384 


0.1394 


0.1375 


0.1434 


0.1424 


0.1315 


0.1345 


0.1355 


0.1345 


0.0000 





The matrix in the lower left was the uncorrected pairwise p-distance among the ND2 sequences; the matrix in the upper right was the 
uncorrected pairwise p-distance amont the COI sequences. 



reconstruction and Bayesian inference were 
listed in table 3. 

The relationships within the ornatipes group 
were not stable revealed by different tree- 
building methods as the low supports for the 
basal nodes (Figures 3, 4 and 5), but it was 
surprising that S. mengla was debarred from 




the ornatipes group in all trees, and it was 
placed at the most basal clade of subgenus 
Steganina receiving great supports (MP BP, or 
bootstrap percentages of the MP analysis = 
100; ML BP = 100; PP or posterior 
probability of the 2-/6-partition Bayesian 
inferring = 1.00/1.00). The remaining species 
of the ornatipes group were recovered as a 

S. angustigena-YH 



1.00/1.00 

00/i.oori.s. angustigena-QX 
S. angnstigena-\\H 
S. xipengi sp.nov 
5. pilosella 
S. lingnonensis 
S. albiventralis 
5I rS. zhaofengi-YH\ 
I r/iao/'eng;-YN3 
S. zhaofengi-\U2 
S. nulliseta 



ornatipes 
group 



S. nigrotimbata 



nigrolimhata 
group 



5. xiaoleiae 
S. ctenaria ~~ shlrozui group 

S. undulata undulata group 



subgenus Steganina 



S. mengla-YNi 



-S. qnadvala 
.*>. emeiensis 



»S. mengla-Ym 1 
J subgenus Slegana 



0.1 substitution site 



Figure 3. Bayesian tree of the ornatipes group deduced from the ND2 and COI sequences. The numbers above the branches 
show the posterior probabilities of the corresponding node in the Bayesian inference. High quality figures are available online. 
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monophyletic group with robust supports in 
all trees (MP BP = 100; ML BP = 100; PP = 
1.00/1.00 in the 2-/6-partition Bayesian 
analyses, respectively). The nigrolimbata 
group appeared to be the closest relative to 
this monophyletic group with well supports 
(MP BP = 100; MLBP = 100; PP = 1.00/1.00) 
(Figures 3, 4 and 5). The Bayesian analysis 
yielded a general topology (Figure 3), which 
was mostly congruent with the result of the 
ML analysis (Figure 4). The monophyletic 
group diverged into two branches. One consist 
of S. zhaofengi triple and S. nulliseta, and the 
other further diverged into S. albiventralls, S. 
lingnanensis, S. pilosella, S. xipengi and S. 
angustigenai triple orderly in the Bayesian 
tree, whereas S. pilosella diverged prior to S. 
lingnanensis (ML BP = 41), leaving S. 
lingnanensis and S. xipngi as sister group in 
the ML reconstruction, but with a low support 
(ML BP = 32). The MP tree (Figure 5) 
differed from the ML and Bayesian tree at 
several points. It suggested a very basal 
position for S. lingnanensis in the ornatipes 



group and S. xipengi clustered with S. 
pilosella which was consistent with the 
Bayesian analysis. The Yunnan (-YN), 
Guangxi (-GX) and Hainan (-FIN) samples of 
S. angustigena clustered together with well 
support (MP BP = 100; ML BP = 100; PP = 
1.00/1.00), and so did in the YN1, YN2 and 
YN3 samples of S. zhaofengi (Figures 3, 4 and 
5). 

To determine the relative contributions of the 
two data partition to the combined analysis 
tree, partition Bremer supports were 
calculated and given in table 5. Support for 
combined analysis phylogeny came from ND2 
was a little bit more than that from COL 
Nodes 1, 3, 6, 7 and 12 showed a mixture of 
positive and negative PBS scores. 

Discussion 

The phylogenetic trees showed that the 

ornatipes group clearly appeared to be 

paraphyletic. To eliminate the effect of 



Kill 
100 



S. angnstigena-YN 
S. angusligena-GX 
S. angusligena-HN 
— S. xipengi sp.nov. 
S. lingnanensis 



- S. pilosella 



• S. albivemratis 
S. zhaofengi-YN I 
$. zhaofengi-YHi 



'S. zhaof'cngi-YN2 
^— S. nulliseta 



orttalipe.' 
group 



nigrotimbata 
group 



- S. nigrolimhata 

S. xiaoleiae 

- xkiruzui group 

- S. undulata undulata group 



S. mengla-YK] 



s. neng/a-YN2 



subgenus Stegattine 



- 5, quadrate 

S. emeiensh 



] 



iubyenus Sceganina 



(I 1 siihstiiuix 



Figure 4. ML tree of the ornatipes group deduced from the ND2 and CO/ sequences. The numbers above the branches show 
the bootstrap percentages (BPs) of the corresponding node in the ML analysis (-InL = 7925.6953 I). High quality figures are 
^available online. 
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ornatipes group 



/( igro lim hula group 

| shirozut group 
| undulata group 



subgenus Stegana 



subgenus Steganina 



Figure 5. Strict consensus tree of two most parsimonious trees of the ornatipes group deduced form the ND2 and COI 
combined sequences. The numbers above the branches show the bootstrap percentages (BPs) of the corresponding node in 
the MP analysis [tree length = 1225, consistency index (CI) = 0.6008, retention index (Rl) = 0.6362]. High quality figures are 
available online. 



individual difference, another sample of S. 
mengla was included in the analysis, but the 
situation did not change. In morphological 
viewpoints, the S. mengla holds the same 
diagnostic characters of the onartipes group, 
which contradict with our molecular 
phylogeny as it formed a separate branch in 
the phylogenetic tree. The amount of genetic 
divergences between S. mengla and other 
species within the ornatipes group were high 
and overlapped to some extent with the 
divergence between species groups. Although 
speculative, the morphological convergence 
should be the reason for this situation. The 
convergent morphological evolution seems to 
be common in the subfamily Steganinae 
(Otranto et al. 2008), which is similar to the 
suggestion made in this research concerning 
convergent morphological evolution in S. 
mengla. Considering the closer relationship of 
S. mengla with the outgroup S. emiensis 
respect to the other species showed in the 



phylogenetic tree, it is possible that S. mengla 
is the interim species of the divergent between 
subgenus Steganina and subgenus Stegana. Of 
course, this hypothesis should be proved with 
analysis of suitable species of both the 
subgenus Steganina and subgenus Stegana. 
Except the S. mengla, the branch consist of 
the rest species of the ornatipes group showed 
the closer relationship with the nigrolimbata 



Table 5. Partition Bremer support for all nodes in the 
MP tree (Fig. 5). 



Node 


Bremer support 


Individual gene partitions 


ND2 


CO/ 


1 


8 


8.50 


-0.50 


2 


22 


15.00 


7.00 


3 


1 


1.50 


-0.50 


4 


1 


1.00 


0.00 


5 


40 


23.50 


16.50 


6 


4 


-1.90 


5.90 


7 


1 


-4.17 


5.17 


8 


13 


7.00 


6.00 


9 


19 


13.50 


5.50 


10 


17 


13.50 


3.50 


1 1 


9 


3.83 


5.17 


12 


4 


-2.50 


6.50 


13 


70 


43.17 


26.83 


Total PBS 


209 


121.93 


87.07 
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group than other species groups of subgenus 
Steganina was consistent with the 
morphological affinity in the two groups (Cao 
and Chen 2008). 

In general, the NADH dehydrogenase subunit 
genes are rapidly evolving, but the 
cytochrome oxidase subunit is more slowly 
evolving (Simon 1994). It was supposed that 
the ND2 gene was better than the COI gene 
suited for species-level analysis, but the PBS 
analysis indicating that the contribution to the 
MP reconstruction in this research of the ND2 
gene was nearly the same as the COI gene. 
Our PBS analysis had implications for the 
conflicts of the two genes at some nodes (e.g., 
nodes 1, 3, 6, 7 and 12), suggesting that the 
two partitions data (ND2 and COT) may be 
favoring an alternative tree topology. The 
relationships of these nodes should be viewed 
cautiously. 

The genetic distances between Yunnan, 
Guangxi and Hainan samples of S. 
angustigena [p-distance of ND2 = 0.0129 (- 
GX vs. -HN), 0.0120 (-HN vs. -YN), 0.001 (- 
GX vs. -YN); p-distance of COI = 0.0044 (- 
GX vs. -HN), 0.0044 (-HN vs. -YN), 0.0000 
(-GX vs. -YN)] were among the mean 
intraspecific variability of Meier et al. 2008 
(1.3 ± 1.6%) for Diptera. In addition, no 
essential morphological character was found 
to distinguish the specimens of these three 
samples, indicating that they should be taken 
as conspecific ones. It was the same as the 
case of the YN1, YN2 and YN3 samples of S. 
zhaofengi. The genetic data [p-distance of 
ND2 = 0.0050 (-YN1 vs. -YN2), 0.0040 (- 
YN2 vs. -YN3), 0.0010 (-YN1 vs. -YN3); p- 
distance of COI = 0.0000 (-YN1 vs. -YN2), 
0.0059 (-YN2 vs. -YN3), 0.0059 (-YN1 vs. - 
YN3)] also indicated the conspecific status of 
the three samples of S. zhaofengi. 



Lu et al. 

Some relationships within the ornatipes group 
were not well resolved, especially the 
alternative placement of S. lingnanensis, S. 
xipengi and 5*. pilosella. Therefore, it may be 
worthy to increase either the genetic markers 
(such as nuclear markers) or the number of 
samples in the future phylogenetic analysis of 
the ornatipes group. 
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